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Abstract 

The aim of this note is to systematize our knowledge about identical configurations of ASM. 

1 Introduction. 

Abelian sandpile model (ASM) was introduced by Bak, Tang and Wiessenfeld in their work [2T] describing formation 
of avalanches. In most general formulation the model can be defined as following automata. Let Q — (Vg, Eg) denote 
a finite graph. For any vertex v £ Vg denote by N(v) the set of all adjacent vertices N(u) = {Vj £ Vg | (v,vj) £ Eg} 
and by deg(w) = |N(w)| the degree of v. 

Fix some positive integer parameter k. Integral- valued function rj : Vg i-> {x + N} is called a configuration on graph 
Q with potential x. 

Sandpile transformation S acts on the space of configurations Conf(Cf) by two steps: 

1. Increase value of tj(vq) ^ t)(vq) + 1 for randomly chosen vertex Vq £ Vg. 

2. If an updated value of rj at some vertex v' exceeds its critical value x + d(v') topple rj at v' i.e. 

• j](v') i-» r)(v') - d(v r ) 

• rj(v) i — ^ r)(v) + 1 for all v £ N(u') 

Such relaxation process may be written in the form 

topple(r/) =r\ — AI(ry(v) > x + d(v)) 

It is natural to set number d(v) to be equal deg(v) so that total norm of the configuration will not change during 
the toppling procedure. However in this case relaxation process described above will never stop for configuration 
i](v) = (x + deg(w)) + S Vo and thus Sandpile transformation will be ill-defined, natural way to avoid this is to define 
a set of boundary vertices dQ = {fbound} were toppling will decrease the configuration rj{v^ oun d) by some number 
^(vbound) > deg(«bound) an d thus total weight of the configuration = viv) wm dissipate through dQ. 

vEVg 

Original situation considered in [21j provides highly illustrative example. Let Q be a bounded subset of two-dimensional 
lattice Z 2 . Every internal node has exactly four neighbours and its degree is also 4. On the other hand any boundary 
node has strictly less then four neighbours. For Conf = A Vg toppling of the node always decrease the value of 
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Figure 1: a) Boundary nodes has less then four neighbours, b) Burning test 



configuration by 4 and so every boundary node dissipate the total weight of configuration each time toppling process 
goes through it (see Fig [I]). 

Theorem 1 ((see [17])). Sandpile transformation is well-defined for each finite graph Q with boundary: S^r] depends 
only on initial configuration r\ and vertex i>£ and does not depend on the sequence in which toppling procedure were 
done. 

Transformation S being very non-local in the sense of Hausdorf metric \<q — rf\u = \{v £ Vg : rj(v) =^ v'( v )}\ on the 
space of configurations, defines thanks to theorem [l] Markov process on Conf with very remarkable properties. 
High interest to this Markov process was caused by the critical behaviour of the distribution of the quantities |<Sr? — r)\<n 
and \Srj\i — \r)\i. Set of recurrent states for the process S was thus very intensively studied over the past decades. In 
this section we state some theorems which can be found in [IB], [H] , [H] , [H] [S] , [T7] and references therein describing 
the structure of this set. 

Definition 1. Let r\ be a configuration on Q with spin k. 

• Vertex v £ Vg is called O-erasable for configuration r\ if n(v) ^ jt + deg(v). Set of all O-erasable vertices for 
configuration n is denoted by £a(r]) 

• Vertex v € Vg is called j -erasable for configuration n if 

n(x) ^ k + deg(w| ) 

Vg\ U £ k (v) 

Set of all j -erasable vertices for configuration r) is denoted by Ej (n) 

N 

• Configuration n is called erasable if there exists such N that Vg = |J £j (rj) . Set of all erasable configurations is 

j=0 

denoted by £(x,Q) 

Theorem 2 (see |17p. Set &{k,Q) coincides with the set of recurrent configurations of the Sandpile process. 
From definition [l] one can easily notice that 

Remark 1 (Monotonicity 1). If n € E(x,Q) and rf > r\ then rj £ £(x,Q). 

Definition 2. We shall say that graph Q\ is embedded in £/ 2 in the sense that Vg C Vg, Eg C Eg and for each v £ Vg 1 : 
deg 2 (w |y s = deg 2 (v). 
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Using again the classical graph for Sandpile model we shall say that subset Oi C Z 2 is embedded into C Z 2 if 
Q.1 C Q, 2 . 

From definition it immediately follows that 

Remark 2 (Monotonicity 2). For rj 6 £(x, G2) the restriction r)\g 1 £ £(x,(/i). 
One of the most remarkable properties of the set E(x,Q) is presented in the next theorem 
Theorem 3 (see [18). Set 8.(k,Q) is bijective to the set of all spanning trees on Q. 
There exists natural bijection £(x, Q) H> £(x + 1, Q). Namely 

?je£(x,e)<^77 + Te£(x+l,e) (1) 

where 1 denotes a function 1 : v € Vg 1— > 1. Thus one can consider £(x, (?) for one value of k. Unfortunately, bijection 
|l| does not hold algebraic structure of the set £,(x,Q) which was noticed in fundamental paper [TT] . 

Theorem 4 (see [11]). Set £(x, Q) with the operation 77, 77' — > r\r{ := topple^ + ?/) is isomorphic to the set J 7 / a of 
equivalence classes of functions on Q up to the image of the Laplace operator. Such a factor space has a structure of 
Abelian group. 

In this paper we are mainly interested in the identical configuration, i.e. erasable configuration which belongs to the 
class of equivalence of {0}. In the section [2] we present some theoretical results concerning identical configurations on 
graphs. Section [5] will be dedicated to experimental results. In section [3] we describe identical configuration for the 
Sierpinskii graph. At last in section [4] we provide a proof of an upper bound of |<St7 — rj\<H on 1? and pose some open 
questions. 

Some results on critical behaviour of Abelian Sandpile Model (ASM) on Z 2 can be found in: [5S], [3], [TS], [3T], [17] . 
|13) and references therein. 

Connection between ASM and similar models on Z 2 is observed in: [10], [5], [19], [26].[2], 

Neutral configurations of ASM on Z 2 were addressed by Creutz in [S] and were studied in [23], [53J, pQ. 

ASM on other graphs such as Sierpinski graph and other self-similar fractal structures were studied in [5], [IB], [7], 

m m- G3- m- m- 

Acknowledgements. Author is deeply thankfull to E.I. Dinaburg and A.N. Rybko for fruitfull discussions. 

2 Some theory. 

Theorem [4] leads to the definition 

Definition 3. Define Green function for two erasable configurations as follows 

Theorem 5. (see 111]) \v) equals the number of topplings occurred at v in the process of relaxation of the 

element rj + 77' 1— > 7777'. 

Proof. Compute the number of incoming and outcoming particles at vertex v in the relaxation process. Total income 
has the form 

G M) (v) -deg(7j)G ( "'"' ) (7j). 

□ 
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Now we are able to notice some properties of the function G^ 1 ' 7 ' \ 

Proposition 1 (Monotonicity I). Ifr}>n' then for any h e £(x,£) it follows that G (?? ' /l) ^ G (?) ' ' h) . 

Proof. Use Theorem [5| If 77 > 77' then relaxation of h + rj can be considered as two consecutive relaxations thanks to 
Abelian property. 

topple(/i + 77) = topple((?7 — 77') + topple(/i + 77 )) 

□ 

Remark 3. Simple computation yields 

AG [r >- h) =rj + h-r)h = rj' + h+ (r) - rj') - r)h = 
= rf + h - r]h + (77 - 77') - (77/1 - /y'/i) = AG iv '' h) + (77 - 77') - (77/1 - 7//1) 

£/l7iS 

r/h-r/'h^r)- 77' (2) 

/or a?7,7/ 77, 77', ft, G £(>f, 5). 
Proposition 2. for otit/ 77, 77' and h 

Q{r),h) _ Q(n'M) _ Q{ri,ri'h) _ ^(jjM') 

Proof. Goes from the definition 

rrn'h = 77 + n'h - AG^^'V =v + r} ' + h- A{G^' h) - G ( "'"'' l) ) = 

fh + 77' + A(G ( '" M ' ) - G^'V - G ( "'"'' l) ) = 77/177' + A(G ( "' i ' ,) ' ) + G ( "'' i} - G (r, ' ,/l) - G (rM '' l) ) 
Thanks to Dirichlet boundary conditions the only harmonic function is identically zero. Which yields the result. 

□ 

Proposition 3. For any 77, 77' 

£ AG^'\v) = J2 deg(v bound )G^'\v bound ) 

VEVg Abound 

Proof. Compute the particles which topples out of the boundary. □ 
Definition 4. Unique configuration Id G £(>r, £/) such that for any other configuration rj € £,{x,Q) 

\Atj = rj 

is called identical configuration. 

Existence and uniqueness of such configuration is granted by Theorem [4] 
Denote by 77* maximal configuration. 

rj* (v) = x + deg(u) 

Theorem 6. 

min max = G (ld/ "* } =: G* 

v v' 



4 



Proof. Proof goes in two steps. First, by theorem jlj max ) = G ,hv . 
At second, since rj* ^ h for any h £ £(x, (?) then by theorem [l]we get 

Since Id G {0} it means that Id = AG^ for some G* £ Zg. Thus G^ 1 ' 1 ^ = G„ for any n € £.(x,Q). In particular, 

□ 

Proposition 4 (Monotonicity II). I/t/i is embedded in Qi then 

Proof. Let Id' denote the restriction of identical configuration Id £ £(^,£2) on £ 2 to the set Vg. Then G^ v ,ld ) ^ 
G^ \ v i. Denote 

h = n* - Id'rf 

Then Gr>*> h = G r '*< ld ' and G''*'' 1 > G^ 1 since obviously n* + h £ {77*}. □ 

3 Identity on Sierpinski carpet. 

As it was mentioned in the introduction, there is a natural bijection ([TJ between two sets of erasable configurations 
with different spins. Unfortunately, in general 1 doesn't belong to {0} and so bijection isn't isomorphic. Thus 
question about identical configuration is the question about the orbit of function 1. 

Proposition 5. For any n £ N there exists such x that Id^ £ {(1)™}. 

Proof. Since every object under consideration is finite there exists a cycle in the sequence {(l)' c }fe e N- Thus the set 
{(l) fc }fceN forms a subgroup in the set of all erasable configurations E(x). For instance there exists such n that 
/ £ {!"} n £(x) = {IT 1 }. Thus applying bijection one gets / + T £ {0} n £(x + 1). □ 

We present the following table illustrating, the fact, that such orbit can be sufficiently large. Consider f2 containing 
only three consequent cells. 

' x I y 



It follows from symmetry that for any x \&{x) = ld(z) so one can get. 
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G(x) 


G(y) 
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5g + l 


6g + l 


14g + 2 


1 





5g+ 1 


6g + l 


14g + 3 
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5g + 2 


6g + 2 


14 9 + 4 


2 





5g + 2 


6g + 2 


14g + 5 
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5g + 2 


6g + 3 


14 9 + 6 
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5g + 3 


6g + 3 


14 g + 7 
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5g + 3 


6g + 4 


14g + 8 
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5g + 3 


6g + 4 


14 g + 9 
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5g + 4 


6g + 5 


14<? + 10 


1 


2 


5g + 4 


6g + 5 


14g + ll 


3 


3 


5<? + 5 


6g + 6 


14g + 12 


2 


2 


5g + 5 


6g + 6 


14g + 13 


1 


1 


5g + 5 


6g + 6 


14g 


3 


2 


5q+ 1 


6g + l 



Remark 4. There are cases with two possible configurations in the table. For x — 14g + 11 one get |0|1|0| which is 
unerasable and |3|3|3| which is erasable. Similarly, for x — \Aq one get unerasable |0|0|0| and erasable |3|2|3|. 



Question about orbit of particular element of the group is very interesting and can be addressed to the future research. 
We shall not cover it in this survey (see section [4]). 

Thus situations when 1 G {0} arc somehow exceptional since in that case question about identical configuration makes 
sense. 

In this section we shall consider another well-known regular graph of order 4 - Sierpinski carpet. 
The only argument to consider such a fractal here is the following 

Theorem 7. On N-th Sierpinski carpet identical configuration has the form 

ld(x) = x + 3 for x = 2n+l 
However for n — In identical configuration does not equal to constant (see Fig [3]). 

We hope that there exists an elegant proof of this fact different from ours which is just constructible. We shall prove 
the existence of the function corresponding to the identical configuration. 

Since all sites in Id has the same value 4 one can easily compute the value of G at the boundary. While toppling sum of 
two identical configurations energy can dissipate only through the boundary points with the rate equal 4— deg(i>bound) = 
4 — 2. From symmetry it follows that energy will dissipate from all three vertices equally. Thus one should only calculate 
the number of particles in identical configuration, which is an easy task. Thus 



G(wbound) = 



i 1 HlcJIU = i ^ 

3 4-2" " 3 4- 



4-3-1 



3 W+1 - 1 
2 



= 3 W+1 



(3) 



Secondly we shall prove reduction lemma 



4 4 

Lemma 1. If AG(v a ) = H + £ G(a i ) where a j e S n then 4G(z ) = 5ff + E G ( b i) with b j € 

i=i j=i 



Proof. Proof consists in straight calculation: 
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Figure 2: Identical configuration on Sierpinski graph for k 



Figure 3: Reduction lemma 



4GM = H + G(ai) - 
4G(a x ) =H + G(v ) - 
4G(a 2 ) = Ho + G(v ) - 
4G(a 3 ) = H Q + G(v ) - 
4G(a 4 ) = H Q + G(v ) - 
4G(a 5 ) = H + G( ai ) 
4G(a 6 ) = H + G(a 2 ) 

Sum over <[5j)~([8]) and introduce Q 



G(a 2 ) - 

- G(a 4 ) - 
■ G(a 3 ) - 

- G(a 2 ) - 
-G(aO- 

- G(a 4 ) 



G(a 3 )- 
■G(a 5 ) 
■G(a 6 ) 
■G(a 6 ) 
■G(a 5 ) 



G(a 4 ) 
■G(6 X ) 

■ G(6 2 ) 
G(6 3 ) 

■ G(6 4 ) 
- G(6 4 ) 



G(a 3 ) + G(b 2 ) + G(b 3 ) 



Sum separately ^ and ( 10 ) 



2^G(a,) = 5i/ + E G ^) 
.7=1 i=i 



2(G(a 5 ) + G(a 6 )) 



4(G(a 5 ) + G(a 6 )) - 2ff + ^ + £ Gfo) 



Multiply (11) by 2 and introduce (12) 



4 

4^G( % ) - lOifo +2^G(&,-) + 2H + J2 + ]T Gfo, 
i=i i=i j=i j=i 



or 



4 4 



£G(a,)=4ff Q + ^G(&,) (13) 
j=i i=i 



Introducing ( 13 ) into Q we obtain 
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From reduction lemma it follows that 



2 

Denote by H n := 5 n H and by M n := G(v N- n ) — G(a) then 



4Gf(«o) = 5F +X) G ( & i) ( 14 ) 

3=1 

□ 



G(«o) = G(a) + V# (15) 



2 

Lemma 2. 



Afjv = \h n (16) 



3 
5 



From lemma [2] it follows 

o oN—n i 

M n = -^— 2 H n = (4 • 3 N - n - 2)5" (17) 

Finally to reconstruct the function G on the whole graph we use the following lemma. 

Lemma 3. For given values G{a), G(b) and G(c) for a, b, c G §™ value at the point v a € §™ +1 which lies on the side 
opposite to a is equal to 

G(v a ) = l(G(a) + 2G(b) + 2G(c) + 3H n ) 
5 

Proof. Using the same trick from lemma [l] we can write 

G(v a ) + G(v b ) + G(v c ) = G(a) + G{b) + G{c) + 2H n 

and so 

AG{v a ) = G{b) + G(c) + G(v b ) + G(v c ) + H n = G{a) + 2G(b) + 2G(c) + 3H„ - G{v a ) 

□ 
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4 One particular result. 

Here we shall consider one particular case of and provide one locality result which can be useful for construction of 
limiting dynamics. 

Denote by V R the diamond of radius R 

V R = {x | |x| Ma n := \xi\ + \x 2 \ «S R} 

Theorem 8 (Locality). For any fixed r and sufficiently large R for any r\ € £,{x,'D R ) such that supp(?7 cr — /} C T> r 
and for any point x E T> r 

supp{(6 x r) cr ) - 8 x rj)} C V r+3 U {ati = 0} U {x 2 = 0} 
4.1 Proof of theorem [8l 

First we deduce the statement of the theorem from some pure constructive proposition and after that we shall present 
proofs of that propositions. 

Proposition 6. For any f from the statement of the theorem and for any point x there exists such configuration f 
that S x f = 5 f. 

So we can consider the most general case x = 0. Denote by / 

x + 2d-3, |x|i =r + 2, \x lX2 \ 7^0 
f(x) = { x+2d-2 7 \x\i =r + 2, \x x x 2 \ =0 
f(x), else 

Proposition 7. / E E(V R ) f E E(V R ). 

Proof. Since all of the points in V R \ T> r+2 are obviously erasable, it is sufficient to show that / G £(x, T> r+2 ). Wc 
will show even stronger result that every point in the belt T> r+2 \ T> r is erasable independently of the configuration / 
(since it differs from 3? only in T> r ). Then since / E £(x, T> r ) the statement will be proven. 

Points (±(r + 2), 0), (0, ±(r + 2)) are 0-erasable by definition. Points (±(r + 1), 0), (0, ±(r + 1)) are then 1-erasable, 
since their value is 2d — 1 + k. Points (±(r + 1), ±1), (±1, ±(r + 1)) are 2-erasable and so (±r, ±1), (±1, ±r) are 3 - 
erasable and so on. 

In general points {xx,x 2 ), \x\\ = r + 2 are (2min(|a;i|, |a;2|))-erasable and their inner neighbours (xi — l,x 2 ) and 
(xi,x 2 — 1) are (2min(|a;i|, j^l) — l)-erasable and (2min(|a;i|, \x 2 \) + l)-erasable consequently. □ 

Thus one can write 

/ = £(±(r+2),0)<*(0,±(r+2)) ^/ (18) 

\«\l = r + 2 



and all operations 5 X in (18) commute so for Sof we get from ( |18[ ) 

3o.f = <5(±(r+2),0)^(0,±(r+2)) J| S y( S of) (19) 



l»ll=r+2 



It is easy to check that 



suppG^ 3<+2d - 1) C V r+2 
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Then from theorem [l] one can conclude that nothing topples out of the region T> r +i for the configurations / for any 
erasable configuration /. 



If for any x : \x\i — r + 2 we have 5of(x) = f(x) then from ( 19 ) Sof(x) < x + 2d — 1 and nothing topples out of the 
T> r+2 so the statement of theorem is satisfied. 

Else there are some points x on the boundary such that 5of(x) > f(x) so after returning excavated particles they 
should topple. In other words for some x one will get 

fiof(x) + <5(±( r+ 2),0) + <5(0,±(r+2)) + 2 ^ S y > X + 2d - 1 

\y\l=r+2 

We shall carefully follow the process of toppling and prove by induction, that any site x of the configuration topples 

not more than R — \x\% + 1 times. 

We shall distinguish two kinds of toppling 

1. Toppling in the domain T>r \ XV+i 

2. Toppling inside XV+i 

Toppling of the first kind. 

Lemma 4. For any connected set M and any point x^ £ dM 

G^ m ^ = l 

Proof. If the set M is connected then for any point x £ M there exists a path from x^ to x. Obviously, if this path 
contains only cells with 2d — 1 + x particles and the starting point of this path topples then each cell should topple. 
So 

G 5 *w*(x) > 1 

The aim is to prove that there is an identity. The proof goes by induction of the area of M. For M containing only 
one cell the statement is obvious. Suppose that the lemma is proven for any connected set M consisting of N cells. 
Consider such set M' = M U {x'} that x^ £ dM' . 

Since M' is connected then x' has not more than 2d neighbours from M. By the induction statement any of this 
neighbours toppled precisely one time and so, since x' ^ M each of them became less or equal than x + 2d — 2 after 
such toppling since they have less than 2d neighbours. 

Cell x' receive not more than 2d particles and so topples one time and distribute 2d particles between its neighbours. 
So any neighbour gets 1 particle and became not more than 2d — 1 + x . □ 

Proposition 8. For any belt M = T>r x \ "Dr 2 such that 2 < R\ — R2 and for any point x : = R\ 

G$*(x) = 1 

and 

8 x x=x-2 ^2 s y - <^(±(iJi),o) - £(o,±(fli)) - 2 ^ S v + *(±(-R2),o) + <*(o,±(fl a )) 

l»ll=«l lvll = «3 

\yiv2\i^ G |yiy2M° 

Proof. For R x — R 2 >2 belt M is a connected belt so, by the lemma ||]G , ^o' 3 *m) = 1. 

Any cell which receive as much particles as much neighbours it have and lose 2d particles. It means that any 
cell which does not belong to the boundary receive and lose equal amount of particles, so it remains 2d — 1 + x. 
Cells at the outer boundary lose 2 + 5(±(Ri),0) + <5(0, ±(i?i)) particles. At last cells on the inner boundary lose 
2 - 8(±(R 2 ), 0) - 6(0, ±(i? 2 )) particles. □ 
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Thus for one toppling of the first kind each cell x <G <9£V+2 gets r} x (T>n \ XV+2) particles. Clearly rj x (T>n \ P r +a) = 

2 + 5{±{ r +2),Q) + £(0,±(H-2))< 

Toppling of the second kind. 

Lemma 5. For any set M define a function 



I(x) = 



Then for any f S £(x, M) 



2d~n x (M), x€dM 
0, e/se 



r»(i,/) _ -1 



Proof Obviously, I(x) = Al(x). □ 

Now calculate the number of particles which any cell (x, y) £ 9XV+1 gets during one toppling of the second kind. It 
gets rj x (V r+ i) particles. In other words any cell at dV r+1 receives 2 — <5(±(r + 1), 0) — (5(0, ±(r + 1)) particles. So the 
number of particles in each point except outer boundary stays unchanged after one step of toppling of the first and 
second kind and so some points at &D r+ 2 remains greater than 4. 

Now we can consider only Dr—i instead of Ur and go another step of induction, and while toppling of the first kind 
any cell on the boundary dV R receive 2 — 6(±(R), 0) — 6(0, ±(i?)) particles. So values on the edges remain unchanged 
and values in the vertices become x + 2d — 3. This circumstance finish the proof of the theorem. 

4.2 Some conjectures and open questions 

• There are several questions arising from ([!]). How does the period of the cycle 1™ depend on the set Vgl How 
does this set distributed in the whole set £? What can we say about asymptotic behaviour of the "dimensionless" 
G 

function — — ? 

x 

One can conjecture that for the sets consisting only of their boundary such an orbit contains all "symmetric" 
configurations. Thus for example periods of this orbit for the few first subsets of two-dimensional lattice are: 





T=4, 






T=3, 









T=14, 



T=2, 



T=ll. 



T=13 



However this is certainly not true for the general case. Thus for the square 3x3 such an orbit contains only 16 
configurations. 

• Another set of questions comes from the correspondence to the spanning trees model. Is it right, that the longest 
tree corresponds to the smallest possible erasable configuration? If it is so, then the minimal weight of erasable 
configuration is asymptotically 2 + x which somehow correlates with the weight of identical configuration. 

• What can one say about the mean level of the cell in Vg over all spanning trees? (Cesaro mean)? 

• How does the structure of graph affects the structure of the orbit 1™ and thus how it is related to the criticality? 
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5 Numerical Experiments 



Numerical experiments provide the evidence of some remarkable properties of the identity configuration on Z 2 . Being 
non-invariant under the change x —¥ k + 1 they still preserve their internal structure and self-similar portraits. We 
claim that such a rigidity is caused by the underlying rigidity of the corresponding functions G ^ and thus study of 
these functions, presented in the section [2] might be of some interest. 




Figure 4: Identity configuration for a = 1 seems to be rotational invariant and self similar. Top raw: identical 
configurations for square, diamond and rectangle for x — \. Yellow points correspond to r\(v) = 4, blue to 3, red - 2 
and black - 1. Bottom raw: identical configurations for octagon, circle and ellipse with eccentricity 0.5. 
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Figure 5: Left to right: Critical configuration for x = 1 on the square 605 x 605 with added identical configuration 
for the square 603 x 603. Identical configurations for the square and diamond for k = 0. 
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